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Shock-heating of stellar envelopes: A possible common mechanism 
at the origin of explosions and eruptions in massive stars 
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ABSTRACT 

Observations of transient phenomena in the Universe reveal a spectrum of mass-ejection prop- 
erties associated with massive stars, covering from Type II/Ib/Ic core-collapse supernovae 
(SNe) to giant eruptions of Luminous Blue Variables (LBV) and optical transients. In this 
work, we hypothesize that a large fraction of these phenomena may have an explosive origin, 
the distinguishing ingredient being the ratio of the prompt energy release E'dep to the en- 
velope binding energy -Ebinding- Using one-dimensional one-group radiation hydrodynamics 
and a set of 10-25 Mq massive-star models, we explore the dynamical response of a stellar 
envelope subject to a strong, sudden, and deeply-rooted energy release. Following energy de- 
position, a shock systematically forms, crosses the progenitor envelope on a day time-scale, 
and breaks-out with a signal of hour-to-days duration and a 10 -10 L0 luminosity. We iden- 
tify three different regimes, corresponding to a transition from dynamic to quasi-static diffu- 
sion transport. For E^ep > ^binding, full envelope ejection results with a SN-like bolometric 
luminosity and kinetic energy, modulations being commensurate to the energy deposited and 
echoing the diversity of Type Il-Plateau SNe. For -Edcp ~ ^binding, partial envelope ejection 
results with a small expansion speed, and a more modest but year-long luminosity plateau, 
reminiscent of LBV eruptions or so-called SN impostors. For Ea cp < -^binding, we obtain a 
"puffed-up" star, secularly relaxing back to thermal equilibrium. In parallel with gravitational 
collapse and Type II SNe, we argue that thermonuclear combustion, for example of as little as 
a few 0.01 Mq of C/O, could power a wide range of explosions/eruptions. Besides massive 
stars close to the Eddington limit and/or critical rotation, 8-12 M Q red-supergiants, which are 
amongst the least bound of all stars, represent attractive candidates for transient phenomena. 
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1 INTRODUCTION 

Stellar explosions, broadly refered to as supernovae (SNe), are un- 
derstood to stem from a sudden release of energy either associated 
with the collapse of the degenerate core of a massive star or from 
the thermonuclear combustion of fresh fuel deep inside the stel- 
lar envelope. Whether one or the other mechanism occurs seems 
to depend on the main-sequence mass of the progenitor star, with 
core collapse occurin g systematically if its value is above ~8 Mq 
dWooslev et"alll2002l. hereafter WHW02). Interestingly, whatever 
the mechanism, the typical kinetic energy of SN ejecta is on the 
order of 10 51 erg, as infer red for example for th e well-studied SN 
1987A (Type II peculiar; iBlinnikov et afll2000h . for SN 1999em 
(Type Il-Plateau, heareafter II-P; Iutrobinll200% . or for the very 
uniform set of events that Type la SNe constitutes dWooslev et al.l 
l2007bh . The cause of this apparent degeneracy in explosion energy 
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is, paradoxically, perhaps not so much tied to the mechanism itself, 
but instead to the rather uniform total envelope binding energy of 
the progenitor stars, on the order of 10 51 erg; anything falling short 
of that leads to a fizzle and no SN display. 

The last decade of observations of such transient phenom- 
ena has shown, however, that the radiative signatures associated 
with SNe (as classified in circulars) are very diverse, from very 
faint to very luminous, from fast-evolving to slow-evolving or fast- 
expanding to slow-expanding. This diversity has been observed lit- 
tle in Type la SNe, with a few peculiar events such as SN 2002ic 
(presence of narrow hydrogen li nes in an otherwise standard Type 
la spectrum; lHamuv et ai]|2003h or SNLS-03D3bb ( possible Type 
la SN from a super-Chandrasekhar white dwarf stan lHowell et al.l 
120061) . In contrast, there has been a rich diversity in explosions as- 
sociated (perhaps erroneously at times) with massive stars and the 
mechanism of core collapse. We have observed 1) Type Ic SNe, 
associated or not with a long-soft 7-ray burst, an d with a stan- 
dard or a very large kinetic energy (SN 1998bw, IWooslev et al.l 



© 2009 RAS 



2 Luc Dessart, Eli Livne, and Roni Waldman 



<D 
XI 
3 

e 

OB 

a 

s 



.O 
< 





-3000-2000-1000 1000 2000 3000 4000 
Days since Maximum 



200 400 600 

Days since Maximum 



Figure 1. Left: Compariso n of absolute-V-band-magnitude light curves for an illustrative and non-exhaustive sample of SNe, SN im postors and/or erupting 



LBVs (violet: SN 1999em,lLeonard et al]|2002al ; lDessart & Hillieill2006l: blue: SN1999gi. lLeonard et alj|2002bl: turq uoise: SN 2005cs. IPastorello etal. 
green: SN1999b r, IPastorello et alJl20QgF light green: n Car. lFrewll2004t yellow: SN 1997bs. IVan Dvk et al.ll2000l : red: NGC2363-V1. brissen etal 
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IPetit et alj|200d) . For each, we adopt the distance and reddening given in the associated references. The time origin is that of maximum recorded brightness. 
Note that all these objects have comparable effective temperatures on the order of 10,000 K, hence comparable bolometric correction, making the comparison 
of their absolut e V-band magn i tude m eaningful. We also show on the left side and in black the absolute visual ma gnitude of the galac tic red-supergiant 
stars studied bv lLevesque et alj 120051 black), as well as the bolometric magnitude of the O star models computed by Martins et al.1 120051 gray: we use the 
bolometric magnitude here since O stars are hot and have large bolometric corrections). Right: Same as left, but now zooming in on the time of maximum 
brightness (the color coding is the same as in the left panel). Notice the stark contrast between SN light curves, associated with shorter/brighter events, 
and erupting massive stars, associated with longer/fainter events. Importantly, notice the overlap between the intrinsic brightness of r] Car and that of the 
low-luminosity Type II-P SN 1999br. In this work, we propose that this diversity of radiative displays may be accomodated by a common, explosive, origin. 



1 19991 : SN 2002ap. iMazzali et al.|[20o3) : 2) a large population of 
Type II-Plateau (II-P) SNe in what seems to be the generic ex- 
plosion of a moderate-mass red-supergiant (RSG) star (e.g. SN 
2005cs. iMaund et alj|20ol ; lUtrobin & Chugaill2008l) ; 3) a grow- 
ing number of low-luminosity SNe that share properties with stan- 
dard Type II-P SNe except for being significantly a nd globally 
less energetic (e.g. SN 1997D, IChugai & UtrobirJl2o"oo| ; SN 1999br, 
IPastorello etai]l2004 OT2006 -1 in M85, whose stat us is ambigu- 
ous. see iKulkarni et al]|2007l ; IPastorello et alj|2007f) . We show a 
sample of V^-band absolute-magnitude light curves of such core- 
collapse SNe in Fig.Q] with representative peak values of —14 to 
— 17 mag and a 100-day plateau duration (best seen in the right 
panel of that figure), hence about 6-10 mag brighter than their pro- 
posed RSG progenitors (shown as black crosses). From this ex- 
panded SN sample, the range of corresponding explosion energies 
has considerably widened, extending above and below the standard 
10 51 erg value. Within the core-collapse SN context, this modula- 
tion is thought to stem from modulations in the energy revival of 
the stalled shock above the nascent proto-neutron star, in turn mod- 
ulated by the stellar-core structure dBurrows et alj|2007bl) . 

The stretching to low explosion energies of potential core- 
collapse SN events is intriguing. For the most energetic explosions 
belonging to points 1 and 2 above, the classification as a SN is un- 
ambiguous. However, some transient events show an ejecta/outflow 
kinetic energy and a peak magnitude that are SN-like, although the 
events did not stem from core collapse (a star is observed at that 
location on post-explosion/eruption images) ; the community call s 
these SN impostors (e.g. SN1997bs; Fig. [Tl IVan Dvk et alj|200oh . 
Conversely, this raises the issue whether low-energy Type II-P SNe 
are associated with core collapse - they might but they need not. 
We illustrate this overlap in radiative properties for a sample of 
such objects in Fig.Q] One such case is the Luminous Blue Variable 
(LBV) r\ Car, whose properties during its 1843 eruption rival those 
of the low-luminosity Type II-P SN1999br. r) Car survived this gi- 



gantic eruption, which shed about 1 of material in what now 
constitutes the homunculus nebula JSmith et alj|2003l) . In contrast 
to core-collapse SNe, such eruptive phenomena in massive stars 
have been associated with the proximity of the star to the Edding- 
ton luminosity I/Edd = incGM/n (n is the mass absorption co- 
efficient). Due to the steep dependence of luminosity to mass (e.g. 
with an exponent of 3.5 for main-sequence objects), this limit is 
easily reached by very massive stars such as r\ Car, or more gen- 
erally massive blue-supergiant stars. In this context, massive stars 
are thought to undergo considerable mass loss when their lumi- 
nosity overcomes the Eddington limit ]] giving rise to a porosity- 
modu lated continuum-driven outflow (Shaviv 2000l lOwocki et all 
1 2004b . Here, this super-Eddington wind constitutes a quasi steady- 
state outflow, and has therefore been thought to be of a fundamen- 
tally different nature from core-collapse SN ejecta. And indeed, 
one refers to a wind for the former and to an ejecta for the later. 
This dichotomy has been exacerbated by the stark contrast in typi- 
cal light curves of eruptive stars (long lived with large brightness) 
and core-collapse SN explosions (short lived with huge brightness). 
In Fig.Q] we show two known eruptive massive stars that highlight 
this contrast. 

However, recent observations may be challenging such a strict 
segregation. First, the recent identification of very fast outflowing 
material ahead of r\ Car's homunculus now suggests that such mate- 
rial was accelerated by a shock, rather than driven in a quasi-steady 
wind, and th us connects the giant outburst to an explosive origin 
dSmithl2008l) . Second, the existence of interacting SNe tells us that 
a massive eruption can occur merely a few year s before explosion. 
For some, e.g. SN 2006gy dSmith et alj|2007l: ISmith & McCravl 
l2007l : IWooslev et al.l2007al) or SN 1994W jDessart et al.l2009l) . the 



1 Note, however, that energy will have to be supplied to the stellar envelope 
to push it over this limit, and in large amounts to explain such a nebula as 
the homunculus. 
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amount is thought to be large enough to decelerate the energetic 
(and necessarily faster-expanding) subsequent ejection. This very 
strict timing of merely a few years, which is orders of magnitude 
smaller than evolutionary or transport time-scales, suggests a con- 
nection between the mechanisms at the origin of the two ejections. 
For SN2006gy, Woosley et al. propose recurrent pair-instability 
pulsations, a mechanism germane to super-massive stars and there- 
fore extremely rare. For lower mass massive stars, this short delay 
of a few years seems to exclude a very-long, secular, evolution for 
the production of the first ejection since this would have no natural 
timing to the comparatively instantaneous event of core collapse. 
Motivated by these recent observations, we explore in this paper 
whether this diversity of events could be reproduced by a unique 
and deeply-rooted mechanism, associated with the sudden energy 
release above the stellar core and the subsequent shock heating of 
the progenitor envelope. This means would communicate a large 
energy to the stellar envelope on a shock crossing time-scale of 
days rather than on a very long-diffusion time-scale of thousands 
of years or more. Although different in their origin, this energy re- 
lease could be a weak analogue of what results in pair-instability 
pu lsations, i.e. a nuclear flash , as identified in the 8-12Mq range 
bv lWeaver & Wooslevl fl979h . 

In this paper, following this shock-heating hypothesis, we use 
ID radiation-hydrodynamics simulations to explore the production 
of explosions/eruptions in stars more massive that ~10M© on the 
main sequence. Rather than focusing on specifi c models, like thos e 
potentially associated with failed supernovae dFrver et al] l2009h . 
we parameterize the problem through a simple energy deposition, 
taking place with a given magnitude, over a given time, and at a 
given depth in a set of pre-SN progenitor star models. We do not 
aim at reproducing any specific observation but, through a system- 
atic ap proach, try to identify important trends, in a spirit similar to 
that of iFalk & Arnettl Jl977t) . However, we depart from these au- 
thors by studying "non-standard" explosions. In practice, we con- 
sider cases where the energy deposited can be both smaller or larger 
than the binding energy of the overlying envelope, but must impera- 
tively be released on a very short time-scale to trigger the formation 
of a shock. Doing so, we identify three regimes, with "standard" 
SN explosions (short-lived transients) at the high energy end, ob- 
jects that we will group in the category SN "impostors" (long-lived 
transients) at intermediate energy, and variable stars at the very low 
energy end. Let us stress here that we do not make the claim that 
all massive-star eruptions, or all transients in general, stem from a 
strong, sudden, and deeply-rooted energy release in their envelope. 
Here, we make this our working hypothesis and investigate how 
much such an explosive scenario can explain observations. We do 
find that this scenario has great potential and should be considered 
as a possibility when examining the origin of massive-star eruptions 
and associated transient phenomena. 

The paper is structured as follows. In !j2] we briefly present 
the stellar evolutionary models of WHW02 that we use as input 
for our ID 1 -group radiation-hydrodynamics simulations. We dis- 
cuss the properties of the progenitor massive stars of WHW02, such 
as density structure and binding energy, that are relevant for the 
present study. We then describe in ^3] the numerical technique and 
setup for our energy deposition study. In ij4] we present the re- 
sults of a sequence of simulations based primarily on the 1 1 Mq 
model of WHW02, discussing the properties of the shocked pro- 
genitor envelope for different values of the strength W. It . the 
depth W,2t and the duration ( S|4.3t of the energy deposition. 
We also discuss in S14.4I the results obtained for more massive 
pre-SN progenitors, ranging from 15 to 25 M© on the main se- 
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Figure 2. Density distribution as a function of Lagrangian mass at the onset 
of collapse for the models of WHW02 evolved at solar metallicity. Notice 
the flattening density distribution above the core for increasing mass. In 
low-mass massive stars, the star is structured as a dense inner region (the 
core), and a tenuous extended H-rich envelope. 

quence. In ij5] we present synthetic sp e ctra computed separately 
with CMFGEN jHillier & MillerJll998l : iDessart & Hilile7l2005bl 
iDessart et al j|2009h for a few models at a representative time af- 
ter shock breakout. In <J(j] we discuss the implications of our results 
for understanding transient phenomena, and in ij7]we summarize 
our conclusions. 



2 INPUT STELLAR EVOLUTIONARY MODELS AND 
PHYSICAL CONTEXT 

2.1 Input Models 

The energy-deposition study presented here starts off from the stel- 
lar evolutionary calculations of WHW02, who provide a set of 
pre-SN massive stars, objects evolved all the way from the main 
sequence until the onset of core collapse. We refer the reader to 
WHW02 for a discussion of the physics included in such models. 
These physically-consistent models give envelope structures that 
obey the equations and principles of stellar structure and evolu- 
tion. The exact details of these models are not relevant since we 
aim at developing a general, qualitative, understanding of stellar 
envelope behaviour in response to shock-heating - we do not aim 
at connecting a specific observation to a specific input. Here, we 
focus on 10-40 Mq progenitor stars evolved at solar metallicity 
and in the absence of rotation. These calculations include the ef- 
fect of radiation-driven mass loss through a treatment of the outer- 
boundary mass-flux, so that the final mass in this set reaches a 
maximum of ~15M© for a ~20M© progenitor star. Mass loss 
leads also to a considerable size reduction of the stellar envelope, 
with partial or complete loss of the hydrogen rich layers, producing 
smaller objects at core collapse for main-sequence masses above 
~30M©. In the sequence of pre-SN models of WHW02, surface 
radii vary from 4-10x 10 13 cm in the 10-30 Mg range, decreasing 
to ~ 10 11 cm in the 30-40 M© range. 

This set of models is not exhaustive. All our input stellar- 
evolutionary models have reached the end of a massive-star life and 
all have a degenerate core that collapses. They are thus not exactly 
suited for post-main sequence massive stars that may be burning 
hydrogen/helium in the core/shell. The sampled mass range does 
not stretch to the most massive stars known, which also happen to 
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exhibit variability and eruptions. In particular, it does not include 
massive blue supergiants, which are connected observationally with 
eruptive phenomena. It does not include highly-massive stars close 
to the Eddington limit, nor objects that for some reason would be 
near critical rotation. It does not include any object resulting from 
close binary evolution. These are clear limitations for an associa- 
tion of transient phenomena with specific progenitors, but this is 
not the goal of this study. Observations have likely been biased 
toward the most massive stars undergoing massive ejections, miss- 
ing fainter progenitors undergoing more modest eruptions. Circum- 
venting specific issues linking tightly specific models to specific 
observations, we attempt instead to develop a qualitative under- 
standing. Our study is thus conceptual and explores the behavior 
of a gravitationally-bound envelope subject to shock-heating. As 
we explain below, the present set covers a wide enough range of 
properties to reveal important systematics in this context, which 
can then be applied more generally to other stars, be they of larger 
mass etc. 

The density distribution above the iron core is an important 
property that distinguishes the structure of massive stars. Low-mass 
massive stars show a steep density fall-off above their degenerate 
core, while the more massive stars boast very flat density profiles 
(Fig. |2j recall that at the time shown, the Fe or ONeMg core is de- 
generate and starts to collapse). In the context of core-collapse SN 
explosions, this has been recognized as one of the key aspect of the 
problem, since this density profile directly connects to the mass- 
accretion rate onto the proto- neutron star in the cri tical first sec- 
ond that follows c o re bou nce dBurrows et alj|2007bl). This has led 
iBurrows & Goshvl Jl993h . and more recently lMurphv & Burrowsl 
1 20081) to identify a criterion for explosion through a competition 
between neutrino luminosity and mass accretion rate. Such a crite- 
rion, although tuned by the global hydrostatic configuration of the 
progenitor star, isolates the conditions for a successful explosion to 
those in the direct vicinity of the pre-collapsed core. A more global, 
and more fundamentally meaningful, criterion for a successful ex- 
plosion is whether the energy deposited at the base of the progenitor 
envelope is smaller or greater than its binding energy. This matter 
is rarely addressed explicitly, supposedly because the 10 51 erg ex- 
plosion energy aimed for is thought to be much in excess of the 
binding energy of the envelope, but this may in fact not be the case 
and it is ignoring the key role the binding energy plays in produc- 
ing the diversity of characters observed in core-collapse SNe, and 
potentially in numerous transients. 

2.2 The Binding Energy Barrier 

The variation in the progenitor envelope density structure (Fig. [2} is 
echoed in the corresponding binding energy of the stellar envelope. 
In the left panel of Fig. [3] we show the binding energy of the enve- 
lope exterior to a Lagrangian-mass coordinate M r , as a function of 
M r , for the same set of pre-SN massive star models. To differenti- 
ate mass shells of different compositions, we use symbols to mark 
the inner edge of the hydrogen shell (diamond), of the He-rich shell 
(triangle), and of the O-rich shell (dot). The binding energy shows 
a strong dependency with progenitor main-sequence mass, and we 
can identify two classes of objects. Objects with main-sequence 
mass below ~25 Mq possess a sizable hydrogen envelope, whose 
fraction of the total mass at the onset of collapse grows as we 
go to lower-mass progenitors. It represents at collapse ~80% of 
the total mass in the 1 1 M Q model, but only ~30% in the 25 M 
model (these masses refer to the main-sequence mass). In all these 
cases, whenever present, the binding energy of the hydrogen shell 



is very small, on the order of 10 erg, hence very loosely bound 
to the star. In contrast, for objects with main-sequence masses in 
excess of 25-30 Mq, the hydrogen envelope is shed through mass 
loss during the pre-SN evolution, and the pre-SN star is essentially 
a "compact" star, with an envelope whose binding energy is now 
3-4 orders of magnitude larger. This is perhaps the single-most im- 
portant characteristics of pre-SN massive star envelopes resulting 
from the single-star evolution calculations of WHW02, namely that 
hydrogen-deficient (generally higher-mass originally) stars are very 
tightly-bound objects while hydrogen-rich (generally lower-mass 
originally) stars are very loosely bound. 

Moving deeper into the envelope, i.e. in the He or O shells, 
the same trend persists, with a systematic increase in binding en- 
ergy for each shell edge as we move up in mass from low-mass 
to high-mass massive stars. As shown in the right panel of Fig. [3] 
the binding energy of the entire envelope exterior to the degener- 
ate core varies from <10 49 erg for the 11 Mq model (a RSG star 
with a final mass of ~1O.6M ) up to ~10 51 erg in the 40 Mq (a 
H-deficient Wolf-Rayet star with a final mass of ~8 Mq) model. 

For an explosion to ensue, an energy greater than the binding 
energy must be deposited suddenly at the base of the envelope, the 
excess energy above the binding energy being used to accelerate the 
envelope and turn it into outflow/ejecta. While a very modest en- 
ergy of ~10 49 erg is sufficient to unbind the envelope of a 10 Mq 
RSG star, such a deposition in a Wolf-Rayet star would produce 
nothing remarkable. In the present work, the connection between 
energy deposition and binding energy is key to the understanding 
of the properties of the resulting mass ejections. More generally, 
the low binding energy of massive star envelopes is of primary im- 
portance for the understanding of their stability/variability, as in the 
context of pulsations, for example. 

The origin of the energy deposition that we artificially treat in 
this work may be a sequel of the gravitational collapse of the core, 
in which case it can happen only once, leading either to an explo- 
sion or to a fizzle. Alternative sources in massive stars will likely 
be less energetic, but in objects with low binding energy, may still 
be of relevance to a wide spectrum of astrophysical events. The 
most suited mechanism in the present context would be the prompt 
thermonuclear combustion of a small amount of material. In the 
context of very massive st ars like r\ Car (not directly addressed 
through our progenitor sefh lGuzikl J2005h propose that non-radial 
gravity mode oscillations could lead to mixing into the hydrogen- 
shell burning layer of fresh fuel located directly above it, thereby 
releasing suddenly a large amount of energy to unbind a fraction 
of the above layers. In lower mass massive stars at the end of their 
lives, this combustion could concern already processed material, 
like carbon, oxygen, or silicon, at a location just above the de- 
generate core. Such "nuclear flashes" have indeed been encoun- 
tered in stellar-evolutionar y calculations of 8-12 Mq massive stars 
dWeaver & Woosley|[l979T) . Interestingly, with a (nuclear) binding 
energy of ~1 MeV/nucleon, the combustion of 12 C or 16 to 56 Fe 
liberates E nuc ~ 1.9 X 10 49 (M/0.01M o )erg. Hence, as little as 
a few percent of a solar mass of carbon or oxygen burnt to iron can 
yield an energy in excess of the binding energy of the lowest-mass 
massive stars, and of compara ble magnitude to that of the weakest 
core- collapse SN explosions JPastorello et alj|2004l ; iKitaura et al.l 
l200d) or LB Vs. 
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Figure 3. Left: Same as Fig. [2] but now for the binding energy as a function of the Lagrangian-mass coordinate M r of the env elope exterior to M r : 
ehmdmg(Mr) = f rtn ^(GM r /r - e int )dM r . In our determination of the binding energy, we have included the internal energy jZeldovich & Novikovl 
Il97ll) . Symbols refer to the base of the corresponding shells, with diamond for the H-rich shell, triangles for the He-rich shell, and dots for the O-rich shell. 
Right: Variation of the envelope binding energy at the onset of collapse in the models of WHW02 and shown as a function of the progenitor mass on the main 
sequence. Here, the "envelope" refers to total mass exterior to the inner 1.8 Mq . The colorbar applies to both left and right panels. 



2.3 Energy Transport by Diffusion/Convection versus 
Shock-heating 

In radiative stellar envelopes, energy is transported outward by dif- 
fusion. Any increase in energy release from the core provokes a 
steepening of the temperature profile and an enhanced radiative 
flux, carrying outward this extra energy, often aided by convection 
too. However, radiative diffusion or convection can only be effec- 
tive for small increases in energy, due to the short mean-free path of 
photons and/or the modest sound speed. Associated time-scales for 
such energy transport are therefore long and the means poorly effi- 
cient. In Fig. [4] we show the diffusion time as a function of depth 
in the 10-40 M© model sequence of WHW02, computed by adopt- 
ing a constant mass-absorption coefficient of 0. lcm 2 g _1 (inter- 
mediate between that for the hydrogen-r ich and silicon-rich sh ells) 
but a depth-dependent mean-free-path dMitalas & SillsllT992L this 
is controlled primarily by the twenty-order-of-magnitude variation 
in density between the core and the surface). Resulting diffusion 
times from regions immediately above the core range from 10 4 to 
10 5 years. In contrast, the shock crossing time-scale of the enve- 
lope is ~ R*/{v B hock), which for Wolf-Rayet to RSG progenitor 
stars with radii of 10 11 to 10 14 cm and even modest, but supersonic, 
shock waves with (ushock) ~1000km s _1 , is typically on the order 
of minutes to days. 

In this study, we want to explore what happens when the de- 
position of energy leads to an energy increase that cannot be reme- 
died either by diffusion or convection transport, but instead has to 
lead to shock formation (even for small energy deposition). This 
leads to a completely different regime since in this situation the 
shock can communicate this extra energy to the entire envelope 
on a short time-scale of days: The stellar envelope hence responds 
immediately to the change of conditions at depth, instead of sec- 
ularly evolving to a bigger/smaller or cooler/hotter star. This di- 
chotomy was recently discussed in the context of the helium flash, 
whose associated en ergy was found to be too small to lead to the 
formation of a shock jMocak et al.l2008ll2009h . In our artificial ap- 
proach, the energy release has to occur on a very short time-scale 
to trigger the formation of a shock. How short will depend on the 
progenitor since, as shown in Fig. [5] the depth-dependent internal 
energy (including both radiation and thermal components) varies 
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Figure 4. Same as Fig. [2] but now showing the variation of the diffu- 
sion time. The c ompu tation assumes unequal photon mean-free-path as in 
iMitalas & s7??s] d 19921) , but a constant and representative mass absorption 
coefficient of 0. 1 cm 2 g — 1 . 

considerably between progenitors, increasing at a given M r with 
main-sequence mass. Forming a shock will require a stronger en- 
ergy deposition, a shorter deposition time, and/or a more exterior 
deposition site in higher mass progenitors. 

2.4 Origin of energy deposition 

In the present study on pre-SN massive-star progenitors, we find 
that an energy deposition that is marginally larger than the envelope 
binding energy yields ejecta that are reminiscent of SN impostors, 
but in the present context, are also reminiscent of the mass ejections 
that must take place prior to core collapse in interacting SNe. The 
short delay between the two ejections, of no more than a few years 
to produce a luminous interacting SN, suggests a shock-heating so- 
lution for both ejections. Any other energy-transport means in a 
massive star that could propel mass out of the gravitational poten- 
tial would take place on a time-scale of thousands of years or more. 
The disconnection between the surface and the core thus suggests 
that both events are tied to a phenomenon that happens in the core 
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Figure 5. Same as Fig. [5] but now showing the variation of the internal 
energy (including radiation and thermal parts). 



or in its direct vicinity. We propose thermonuclear flashes associ- 
ated with shell burning in the last stages of massive star evolution 
(pair-instability pulsations are analogous to this scenario but may 
only app ly to stars with a ma in-sequence mass in the range 95- 
130 M : IWooslev et al]2007 et) . These burning stages are very short 
(month to year time-scales) and always occur at the end of the life 
of all massive stars. They offer a natural tuning and a reproducible 
circumstance for the production of interacting SNe, with at least 
one major interaction between the penultimate shell ejection and 
that resulting from the gravitational collapse of the degenerate core. 
Given the relatively large number of interacting SNe compared to 
the strict time requirements to produce them (a few-year delay is no 
more than an instant relative to the lifetime of a star!) suggests that 
indeed the conditions must be met often rather than encountered by 
chance . The stellar-evolutionary calculations of lWeaver & Wooslevl 
dl979h for 8-12 M© massive stars support the occurrence of nu- 
clear flashes in association with Ne and Si core/shell burning, and 
occuring merely a few years prior to core collapse. These simula- 
tions are old and would need to be revisited in a thorough fash- 
ion, with full hydrodynamics, multi-dimensionality, and high reso- 
lution. 

In the context of more massive stars like r\ Car. lGuzikl J2005h 
propose that non-radial gravity mode oscillations above the core 
and near the hydrogen-burning shell could lead to mixing of 
hydrogen-rich material downward into hotter denser layers. Com- 
bustion of this fresh fuel would yield a burst of energy trigger- 
ing mass ejection and a bolometric brightening. In this work, we 
emphasize the dire consequences of any such energy release on 
loosely-bound massive stars, potentially triggering single eruptions 
as evidenced i n transients or multi ple eruptions as evidenced in in- 
teracting SNe dPessart et ai]|2009h . 

More generally, the hydrodynamical inves tigations of the 
late st ag es of nu clear burning in m assive st ars by Bazan & Arnettl 



dl994 Il998h 



iMeakin & Arnett 



Asida& ; Arnettl feOOOh ; I Arnett et alJ ~T2005) 



J2006I) reveal the multi-dimensional and time- 
dependent nature of the problem, issues that have been largely 
overlooked in stellar evolutionary calculations so far. In particular, 
they find that considerable energy can leak out of the convectively- 
burning shells in the form of waves, with a potential to alter the 
otherwise quasi-steady hydrostatic configuration of the stellar 
envelope. Similar investigations for lower-mass massive stars are 
highly needed. 

These works provide a physical motivation for our investi- 



gations. For the simulations presented here, we do not, however, 
build a fully-consistent picture of the problem, but instead start off 
by imposing the modalities of such an energy deposition. In prac- 
tice, we follow a very simple approach, adopting a constant energy- 
deposition rate for a given duration. We thus neglect any feedback 
of the dynamics on the mechanism causing the energy release. 



3 NUMERICAL TECHNIQUE AND COMPUTATIONAL 
SETUP 

The numerical simulations presented in this work were all per- 
formed with Vulcan/ID (V1D), a one dimensional radiation- 
hydrodynamics code that works in the framework of Newtonian 
physics. It incorporates many methods and techniques, using both 
Lagrangian and Eulerian schemes for planar, cylindrical, and spher- 
ical geometries. The Lagrangian hydrodyna mics scheme is based 
on the staggered-grid method described in iRichtmver & Mortonl 
which uses artificial viscosity for shock waves. So me test 
proble ms of an early version of this code are described in iLivnej 
dl993h . 

For the treatment of non-LTE radiative-transfer problems we 
have recently implemented several solvers with different levels of 
approximations for the radiation fluxes. The coarser method is our 
gray flux-limited-diffusion method, which has also a multi-group 
extension. We have also constructed more accurate schemes for the 
radiation fields using moment methods and full angle-dependent 
solvers, whi ch are similar in nature (but not in details) to the scheme 
described in lEnsmanl <ll991. 1994). Those however were not imple- 
mented in this work. The radiation field is coupled to matter in 
a fully implicit fashion, which guarantees stability and large time 
steps. Since the important physics relevant to this study occurs at 
large optical depth, the multi-group capability of the code is not 
used. Opacities/emissivities are interpolated from a set of tables 
prepared with a separate program. It computes the populations for 
a large number of atomic levels (up to 1000 per species), under 
the assumption of LTE, for typically fifteen species (H, He, C, N, 
O, Ne, Na, Mg, Si, S, Ar, Ca, Cr, Fe, and Ni), and for a set of 
compositions representative of the pre-SN envelope of a 15 M© 
main-sequence star (we generate tables for compositions corre- 
sponding to a mean atomic weight between 1.3 and 25 atomic-mass 
units). For each species, up to 18 ionization stages are included. At 
present, we account for scattering opacity (due to electrons), and 
absorptive opacity due to bound-free and free-free processes. Al- 
though lines can be included, we have neglected them in this study. 

The equation of state takes into account excitation/ionization 
to all available levels of the different atomic/ionic species in the 
composition. The distribution function of these l evels is computed 
using a method similar to the one described in iKovetz & Shavivl 
jl994[) . The resulting electron density is then used together with 
the temperature in order to extract the pressure, energy, chemi- 
cal potential and their derivatives respective to the electron density 
and temperature from a table computed in advance by solving the 
Fermi-Dirac integrals. The pressure, energy and entropy of the ions 
are then added as an ideal gas. 

All our simulations use as inputs the massive-star models of 
WHW02, evolved from the main sequence until the onset of core 
collapse. At present, V1D adopts the same grid as the one read-in 
from the chosen input model of WHW02. The resolution is there- 
fore quite low at the progenitor surface, which affects the accuracy 
of the breakout signal. The range of mass encompassed by the V1D 
grid is set by the inner mass cut where energy is deposited, the 
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excised region being shrunk to a point mass at the origin. We do 
not compute the explosive nucleosynthetic yields, and thus do not 
alter the original composition of the massive-star envelopes com- 
puted by WHW02. In particular, there is no production of 56 Ni and 
no associated heating through radioactive decay accounted for in 
our simulations. In the main parameter study presented here, we 
use the 11M model (named si 1.0 by WHW02). We first ex- 
plore the evolution of the progenitor envelope after depositing an 
energy -Edcp between 5 x 10 48 erg and 1.28 x 10 51 erg at a mass cut 
Af cu t = 1.8Mq (uniformly deposited in the range 1.8 to 2.3Mq) 
and for a duration of 10 s ( i]4.1b . We also explore the sensitivity 
of the outcome to the location of the energy deposition, from a 
Lagrangian mass coordinate of 1.8 to 3, 7, and 9Mq ( i]4.2b . and 
to the duration of the energy deposition, from ten seconds, to one 
hour, one day, one week, and one month ( i]4,3l >. Finally, we perform 
a few simulations using more massive progenitors, with 15, 20, and 
25 M© main-sequence mass ( i]4.4t . A summary of the model ini- 
tial parameters as well as quantities characterizing the main results 
is shown are Tables [T] and [2] Note that in the text, when we men- 
tion the masses associated with the WHW02 models, we mean the 
main-sequence masses, as in the 1 1 M© or the 25 Mq models; 
these do not correspond to the star mass at the onset of collapse, 
which is shown for all models in, e.g., Fig.|2] 




12 3 4 

Time after Energy Deposition [d] 



Figure 6. Time evolution of the mass-integrated total (black), thermal (or- 
ange), radiative (blue), kinetic (green), and gravitational (red; its absolute 
value) energies for the model s 1 1 _5. Notice the interplay between radia- 
tive (which dominates the internal energy over the thermal component) and 
kinetic energy. The time of shock breakout is 2.6 d, at which the total inter- 
nal (kinetic) energy is 8.3 X 10 49 erg (6.7 X 10 49 erg). The time origin is the 
start of the simulation. 



4 RESULTS FOR THE 11 M Q SEQUENCE 

4.1 Effect of varying the Energy Deposition Magnitude 

In this section, we present the results for the sequence of simula- 
tions based on the 1 1 Mq model for a range of energy deposition 
magnitudes. A compendium of parameters describing the set up 
and the results is given in Table Q] Note that, initially (prior to en- 
ergy deposition), the envelope (total) internal energy is ~ 10 49 erg, 
and its gravitational binding energy is -2.2 x 10 49 erg in this 1 1 Mq 
model (the model mass at the time of collapse is 10.6 Mq). 

In our simulations, energy deposition leads to a strong in- 
crease in internal energy (temperature, pressure), leading to the 
formation of a shock which moves outward. In all runs, the shock 
crosses the entire envelope and eventually emerges at the surface 
of the progenitor, giving rise to a shock-breakout signal (more pre- 
cisely, a radiative precursor precedes the emergence of the shock, 
but one generally refers to this whole phase as shock breakout). 
The shocked envelope left behind has been turned into a radiation- 
dominated plasma, with an extra energy that varies from a small 
(E dep = 5 x 10 48 erg in the si 1_0 model) to a large value 
CEdep = 1-28 x 10 51 erg in the si 1.8 model). 

After energy deposition, all the stored energy is internal. As 
the shock forms and moves out, that stored internal energy is con- 
verted into kinetic energy behind the shock, and as the shock pro- 
gresses outward, more and more material is shock-heated and the 
total internal energy goes up at the expense of the kinetic en- 
ergy. As time progresses further, the kinetic energy eventually in- 
creases due to radiation work on the ejecta. Figure [6] illustrates 
this evolution until just after shock breakout for the case in which 
E dcp = 1.6 x 10 50 erg (model sll.5). 

As the shock approaches the surface, the envelope properties 
become qualitatively and quantitatively different, depending on the 
adopted -Edcp value. For a large £"dcp, we find that at shock emer- 
gence, the total envelope energy is in rough equipartition between 
kinetic and internal. As time progresses, all the energy is eventually 
converted into kinetic as radiation, nearly entirely trapped at large 



optical depth, does work to accelerate the ejecta to its asymptotic 
velocity. 

In contrast, for small -Ed op , the kinetic energy of the envelope 
at breakout and later is very small. Here, the stored internal energy 
is weakly enhanced after energy deposition (i.e. Ea cp < Ei n t orig- 
inally), and this excess energy is essentially all exhausted in doing 
work against gravity, merely lifting off the envelope from the grav- 
itational potential well. At shock emergence, the envelope kinetic 
energy is thus negligible. 

Hence, while in all cases from low to high energy deposition, 
a shock forms and eventually emerges at the progenitor surface, 
the quantitative differences are very large between cases. Transi- 
tioning from dynamic to quasi-static diffusion regimes of energy 
transport, we identify three situations: 1) The energy deposition is 
much larger than the binding energy and a SN-like explosion en- 
sues; 2) the energy deposition is on the order of the binding energy 
and a low-energy explosion/eruption results; and 3) the energy de- 
position is lower than the binding energy and the excess energy 
merely shifts the star to a new quasi-equilibrium state from which 
it relaxes on a very long time-scale. In the following subsections, 
we describe each of these regimes individually. 

4.1.1 -Edcp > -Ebinding 

This regime applies to our simulations with Ed ep between 
4xl0 49 erg (model sll.3) and 1.28xl0 51 erg (model sll_8). The 
separation between this sample and cases with lower energy de- 
position is drawn at the threshold value of Ed ep for which the 
entire envelope is ejected. In Table Q] we give a census of the 
properties for this set. Although qualitatively similar, these simu- 
lations give rise to 8.8 Mq ejecta with a mass-weighted average 
velocity in the range 500-3500 km s~\ and kinetic energy in the 
range 2.8xl0 49 (E dop = 4 x 10 49 erg) up to 1.276x 10 51 erg 
(Edep = 1-28 x 10 51 erg). For the former, the binding energy takes 
a visible share of the energy deposited so that the asymptotic kinetic 
energy is sizably smaller than -Edcp. For the latter, the energy de- 
posited is overwhelmingly large and represents 99.7% of the final 
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Table 1. Summary of the parameters and key results for the sequence started with the 11 Mq progenitor model (pre-SN mass of 10.6 Mq). 
Note that all models are not run for the same total time (i cn d), owing to numerical problems at very late times as the density and temperature 
drop in the outer ejecta. Numbers in parenthesis refer to powers of ten in the unit shown in the first column of the same row. The effective 
mass on the V1D grid depends on the adopted inner mass cut Af C uf The time origin corresponds to the start of the simulation, when energy 
is deposited for t e d ep seconds. Starting from the top row down, we give the total energy deposited (in 10 5<) erg), the inner mass shell where 
it is deposited and the duration of that deposition (in s), the time at the end of the simulation (in days), the sum of the mass of all shells that 
at the end of the simulation have a velocity larger than the local escape speed •J 2GM r /r, the ejecta mass-weighted average velocity (in 
km s — 1 ), the kinetic energy of the ejecta (in 10 50 erg) at the end of the simulations, the maximum ejecta velocity (in km s — 1 ), the time (in 
days) and the peak luminosity (in Lq) of shock breakout , the average speed of the shock (in km s — 1 ), the kinetic and the internal (thermal plus 
radiation contributions) at the time of breakout (in erg), the time of peak luminosity in the post-breakout plateau (in days) and the corresponding 
luminosity (in Lq), the duration of high brightness for the transient (time during which the luminosity is greater thanl/50th of L pea k, plateau), 
and the time-integrated bolometric luminosity (in 10 49 erg; note that the time interval varies between models). [See text for discussion.] 



Model Name 


si 1.0 


sll.01 


sll.l 


si 1.2 


si 1.3 


sll.4 


sll_5 


si 1.6 


si 1.7 


si 1.8 


E dcp (10 50 erg) 


5.0(-2) 


7.5(-2) 


1.0(-1) 


2.0(-l) 


4.0(-l) 


8.0(-l) 


1.6(0) 


3.2(0) 


6.4(0) 


1.28(1) 


A/cut (Mq) 


1.8 


1.8 


1.8 


1.8 


1.8 


1.8 


1.8 


1.8 


1.8 


1.8 


^edep (s) 


10 


10 


10 


10 


10 


10 


10 


10 


10 


10 


*cnd (d) 


730 


193 


730 


730 


160 


730 


117 


103 


730 


686 


^ejected (Mq) 


0.13 


5.95 


7.75 


8.60 


8.72 


8.79 


8.80 


8.80 


8.80 


8.80 


{v) m (km s' 1 ) 


42 


76 


121 


296 


510 


789 


1165 


1687 


2412 


3432 


(£kin)end(10 50 erg) 


2.30(-5) 


4.00(-3) 


1.33(-2) 


9.25(-2) 


2.79(-l) 


6.74(-l) 


1.47(0) 


3.08(0) 


6.30(0) 


1.276(1) 


(%ax) (km s _1 ) 


60 


410 


560 


1310 


1800 


2740 


4100 


5970 


8740 


12700 


*SBO (d) 


54.7 


32.1 


21.4 


9.6 


5.8 


3.8 


2.6 


1.8 


1.3 


0.9 


^peak.SBO ( L o) 


1.7(5) 


3.3(6) 


5.2(7) 


1.2(9) 


4.9(9) 


1.6(10) 


4.2(10) 


1.1(11) 


2.5(11) 


5.7(11) 


(■"shock) (km s _1 ) 


86 


147 


221 


492 


812 


1238 


1815 


2622 


3717 


5244 


(Skin)sBO (erg) 


2.1(46) 


2.1(47) 


6.7(47) 


4.4(48) 


1.3(49) 


3.1(49) 


6.7(49) 


1.4(50) 


2.8(50) 


5.8(50) 


(^int)sBO (erg) 


5.2(48) 


4.6(48) 


3.8(48) 


7.5(48) 


1.7(49) 


3.9(49) 


8.3(49) 


1.7(50) 


3.5(50) 


7.1(50) 


tpeak, plateau (d) 


208 


198 


243 


115 


97 


82 


73 


64 


55 


49 


^peak, plateau (L©) 


3.6(5) 


1.5(6) 


2.9(6) 


1.2(7) 


2.9(7) 


6.0(7) 


1.1(8) 


2.0(8) 


3.6(8) 


6.2(8) 


A^plateau (d) 


675 


517 


341 


287 


>155 


180 


>115 


>106 


111 


100 


/ L bol dt(10 49 erg) 


6.6(-3) 


5.6(-3) 


2.5(-2) 


6.7(-2) 


1.30-1) 


2.5(-l) 


4.0(-l) 


6.6(-l) 


1.1(0) 


1.8(0) 



kinetic energy. In all cases, the initial prompt deposition of energy 
turns the stellar envelope into a true "ejecta", and this material is 
lost from the star without any noticeable fallback. 

We find that the shock-crossing time through the envelope is 
between 0.9 and 5.8 d (the distance between the deposition site and 
the progenitor surface is 4.7 x 10 13 cm), with a time-averaged shock 
velocity between ~5000 and 800 km s~ 1 (in the same order from 
high to low -Ede P )B This therefore modulates the delay between 
the energy release (which may be a sequel of gravitational collapse 
or otherwise) and the time of shock breakout, which is the first de- 
tectable signal testifying for the violent energy deposition that took 
place at depth. Importantly, this breakout signal is the unambiguous 
signature that a shock wave emerged from the progenitor surface, a 
signature that would directly eliminate steady-state wind solutions 
for the origin of the subsequent outflow/ejecta. 

Owing to the sizable range in shock speeds, the breakout sig- 
nal varies considerably in duration. From a random walk argument, 
it is born when the decreasing envelope optical-depth r at the shock 
depth Ar below the surface eventually makes the diffusion time 
tdiff = rAr/c shorter than the shock-propagation time to the pro- 
genitor surface t s hock = Ar/i> s hock. This occurs at an optical depth 
r ~ c/n a hock, and can be very large for a small shock-propagation 
speed. As shown in Fig.[7](note that the luminosity in this figure is 
normalized to the peak value to better illustrate its time variation), 
we find shock-breakout durations from ~0.5 up to ~2 hr in the se- 
quence (we refer to the intrinsic breakout duration, hence neglect 



2 Note that the shock is always supersonic and may propagate at very dif- 
ferent speeds depending on location and initial energy deposition. 




0.0 0.2 0.4 0.6 0.8 1.0 
Time since Shock Breakout [d] 



Figure 7. Time evolution of the normalized intrinsic bolometric luminosity 
around shock breakout, as computed with V1D, and revealing the hour- 
to-day range in duration of the shock-breakout signal, caused here by the 
change in shock speed (or explosion energy). Light-travel-time delays are 
not accounted for. The variation in the breakout-peak luminosity is shown 
in Fig. [8] A color coding is used to differentiate models. 

light-travel time effects which would broaden the signal seen by a 
distant observer over a duration of at least ~ i?*/c ~1 hr). Sim- 
ilarly, the peak luminosity at breakout varies enormously with the 
energy deposited or the shock strength, as shown in the left panel of 
Fig. [8] where the time axis is the time since energy deposition. To 
a higher energy-deposition strength corresponds stronger shocks, 
shorter shock-crossing times (earlier shock emergence), and both 
shorter-duration and greater peak breakout luminosities. For the 
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range of models discussed in this section, we obtain values between 
4.9x 10 9 up to 5.7 x 10 11 L Q , a contrast of about a hundred which 
is comparable to the corresponding contrast in asymptotic ejecta 
kinetic energy. So far, the systematics of shock breakout have been 
tied to variations in the progenitor properties such as surface ra- 
dius or atmospheric scale height, while obviously the explosion 
energy, if stretched to low enough values, can lengthen consider- 
ably the breakout signal. But importantly, while the delay between 
energy deposition and breakout varies by a factor of a hundred in 
this model set, it remains on a time-scale of days and hence much 
shorter than the diffusion time on the order of thousands of years 
for the corresponding layers (see i]2.3l >. 

These properties at and subsequent to breakout are visible in 
the conditions at the photosphere, which we show in Fig. [9] for 
the velocity Vphot, the radius i? p hot, the temperature T p hot> and 
the Lagrangian mass coordinate M p hot • Because of the infrequent 
writing of model quantities to file, the breakout phase is not well re- 
solved. This matters little for the velocity, the radius, and the mass 
at the photosphere, which all evolve slowly. However, the photo- 
spheric temperature shows a peak whose strength is underestimated 
in the high energy-deposition cases, characterized by short break- 
out durations (note that the low spatial resolution at the progeni- 
tor surface prevents an accurate modeling of this breakout phase 
in any case). Hence, at breakout, we find an approximate range of 
photospheric temperatures between 5xl0 4 to multiples of 10 5 K. 
In scattering-dominated atmospheres, the color temperature of the 
emergent spectral-energy distribution corre sponds roughl y to the 
gas temperature at the thermalization depth jMihalasll 19781. p 149) 
which is up to a factor of two greater than Tpw l lDessart & HillieJ 
l2005ah . For a thermal (blackbody) emitter at T, the peak of the 
spectral-energy distribution (SED) is at ~ 2900/T4 A (where T4 
is T/10 4 K). At breakout, the spectral-energy distribution for the 
present set of models will thus peak in the range 100-600A. It is 
at breakout that a fraction of the material at the progenitor sur- 
face gets accelerated to large velocities. In the sequence of models 
si 1_3 to sll_8, we find increasing maximum ejecta speeds from 
1800 km s _1 up to ~ 13000 km s~ 4 , typically a factor of four 
higher than the mean mass-weighted velocity of the ejecta. 

The general evolution of the long-term, post-breakout, pho- 
tospheric conditions is qualitatively similar for this set of models. 
After breakout and for a few days, the internal energy stored in the 
shock-heated stellar envelope does work to accelerate the ejecta 
to its asymptotic velocity. For larger energy deposition, the expan- 
sion rate is greater and leads to more efficient cooling of the ejecta, 
thereby mitigating the impact of the larger internal energy initially 
provided by the shock. This cooling is quasi-adiabatic since little 
radiation is lost from the photon decoupling layers - it all occurs at 
large optical depth. After a few days, the ejecta velocity increases 
monotonically with radius (homologous expansion is only approx- 
imately attained since the progenitor radius may represent up to 5— 
10% of the current photospheric radius at all times). As the ejecta 
expand, cool, and recombine to a lower ionization state, the photo- 
spheric velocity decreases with time. The photospheric radius first 
increases, reaching values in excess of 10 15 cm for the high energy 
cases, but eventually decreases due to recombination and the dimin- 
ishing optical depth of the ejecta. The location of the photosphere is 
that of full ionization, which for the H-rich composition of this stel- 
lar envelope occurs around 5000 K. Hence, the photospheric tem- 
perature, after its initial peak at breakout, slowly decreases to level 
off at ~5000 K. Ultimately, the photosphere reaches the inner re- 
gion of the ejecta and the "event" enters its nebular phase. This 
transition occurs earlier for higher values of i?dc P , after about 100 



days in model si 1_8 but up to ~200 days in the s 1 1_3 model (since 
the code crashed for that model, this is an estimate based on the no- 
ticeable trend with energy deposition for that quantity; see TableQJ. 

This evolution in photospheric properties parallels the bolo- 
metric evolution of the explosion, which we show in Fig. [8] In all 
cases, the early breakout peak in luminosity is followed by a fading 
(initially from radiative cooling at the shock-heated surface layers) 
followed by a sustained brightness as the ejecta expand (the rate 
of increase in photospheric radius compensates the rate of decrease 
in photospheric temperature): This is the so-called plateau phase 
of Type II-P SNe. Peak luminosities, reached at a later time for a 
lower energy deposition (from 50 to 100 days after breakout), are 
in the range 3 x 10 7 up to 6x 10 8 L©. Similarly, the plateau dura- 
tion lengthens with lower energy deposition, ranging from 100 up 
to about 200 days (recall that this occurs in models si 1_3 to si 1_8 
for the same ejected mass of ~8.8Mq). At the end of our simula- 
tions, we find time-integrated bolometric luminosities in the range 
10 48 -10 49 erg, which represent, in the same order, l/20th down to 
l/70th of the corresponding asymptotic ejecta kinetic energy. 

Hence, at the high energy end, we obtain light curves that are 
reminiscent of Type II- P SNe (Fig. [TV rang ing from the energetic 
events like S N 2006bp llDessart et al1 l2008[), to the more standard 
SN 1999em ^Leonard et alj|2002al) , through to the low luminosity 
SN 1999br jPastorello et alj|2009h . The luminosity contrast of a 
factor of ~20 in this s et of models is on the o rder of that inferred for 
Type II-P SNe (Fig.m lPastorello et al.l2009l) . For moderate to high 
Edcp values, the plateau durations we obtain are on the order of 
those observed for Type II-P SNe, although we find a lengthening 
for lower Edep values that has not been observed definitely. For 
example, the observations during the plateau phase of SN1999br 
are unfortunately truncated after 100 days, at a time when the SN 
was not showing a clear sign of fading. The subsequent observation 
250 days later, when the SN has faded, does not exclude 99br kept 
its plateau brightness for a longer time, perhaps for a total of 150 
or 200 days (see Fig.QJ. 

To summarize, the regime of explosions described in this sec- 
tion corresponds to that of type II-P SNe, with a range of explosion 
energies and bolometric luminosities that are in line with inferences 
from observations. The rather large energy-deposition magnitude 
(i.e. 4 x 10 49 erg) in this regime favors a scenario involving the 
gravitational collapse of the stellar core and, given the resulting 
complete envelope ejection, the formation of a neutron star rem- 
nant. We wish to emphasize that there is no fundamental limita- 
tion for preventing explosions in this 1 1 Mq model for energies as 
low as a few times 10 49 erg, as in model si 1_3. The corresponding 
light curve would be reminiscent of that of SN 1999br, only slightly 
fainter. 



4.1.2 E, 



dep 



-Ebind 



In this section, we continue the sequence toward lower energy 
deposition, entering a regime where its magnitude is on the or- 
der of the binding energy of the envelope, and, incidently, on the 
same order as the total internal energy in the envelope. In contrast 
with results for high Edcp values, only a fraction of the progen- 
itor envelope is ejected in models sll_01 to sll_2, ranging from 
5.95 up to 8.6 Mq. Despite the incomplete success for explosion, 
all models go through a similar sequence as that of powerful ex- 
plosions. A shock crosses the progenitor envelope, heating it and 
propelling it outward. The shock eventually breaks out, after 10- 
30d (time-average shock speeds of 150-500 km s~ 4 ), with an ac- 
companying peak luminosity of 3x 10 6 up to 10 9 L©. The fastest 
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Figure 8. Same as Fig. [7] but this time showing the un-normalized bolometric luminosity. In the left panel, we use a logarithmic scale for the abscissa, 
showing the time since the start of the simulations. In the right panel, we use a linear scale for the abscissa, showing the time since shock breakout. Stronger 
explosions show shorter and stronger breakout signals, followed by a brighter plateau brightness but shorter plateau length. At intermediate values of the 
energy deposition, the plateau-brightness is followed by a ledge, which corresponds to the time when the photosphere recedes into the more tightly-bound 
(and slowly-expanding) helium-rich part of the envelope - no artificial nor physical mixing occurs in our ID simulations, a process that could smear out such 
a feature. Note that for cases in which E^ cp ~ ^binding! me breakout luminosity is barely larger than the following plateau, which we find here can last for 
up to two years in this 1 1 Mg progenitor star. 




1 10 100 200 400 600 
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Figure 9. Top left: Time evolution of the velocity at the photosphere for the sequence of models associated with the 1 1 Mg model of WHW02. Only the 
energy deposited during the initial ten seconds after the start of the simulation differs between models, all initial model characteristics being the same. A color 
coding is used to distinguish the different cases, with total energy deposition encompassing the range 5 X 10 48 (black) up to 1.28 X 10 51 erg (red). Top right: 
Same as top left, but now for the radius at the photosphere. Bottom left: Same as top left, but now for the temperature at the photosphere. Note that the time 
of breakout is poorly sampled (dumps written every 10000 s), so that the temperature shown at breakout is an underestimate. Bottom right: Same as top left, 
but now showing the Lagrangian mass coordinate of the photosphere as a function of time. Note that the mass interior to the inner boundary of the V1D grid 
is 1.8 Mg for this sequence of models, and is not accounted for in this last panel. 
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ejected material is accelerated modestly, from 400 (si 1 _0 1 model) 
up 1300 km s (sll_2 model). At breakout, the kinetic energy of 
the envelope is a diminishing fraction of the total envelope energy, 
which is primarily stored internally. Asymptotically, the total ki- 
netic energy of the material ejected (a fraction only of the total) is in 
the range 4xl0 47 up to 10 49 erg, with mass-weighted average ve- 
locities in the range 70-300 km s _1 . The long-term light curves are 
characterized by peak luminosities in the range of 10 6 — 10 7 Lq and 
very long plateau durations of 1-2 years, yielding low-luminosity 
long-lived transients (Fig. [8}. Hence, these events look somewhat 
like SNe because they attain luminosities/expansion-rates that are 
comparable to those of low-luminosity/energy Type II-P SNe. But 
they also look like giant eruptions as seen in the most massive lu- 
minous stars, characterized by super-Eddington luminosities on the 
order of 10 7 Lq, with ejected shells weighing few solar masses and 
outflowing velocities of a few hundred km s _1 . The Eddington lu- 
minosity for this 10.6 Mq pre-SN star is L-^dd = iircGM/it ~ 
4.07 x 10 Lq (we use a representative mass-absorption coefficient 
k — 0.34 cm 2 g -1 for this hydrogen-rich stellar envelope), and 
hence, the luminosities quoted above in the range 10 6 -10 7 Lq are 
all significantly super-Eddington. By contrast to radiative-driven 
outflows, the super-Eddington nature of the luminosity is only a 
feature in the current context since the ejecta is already unbound 
after shock passage (this does not exclude secondary effects that 
could stem from the high luminosity). 

The outcome of such a weak energy "explosion" can be two 
things. If the origin is a very weak energy release following core 
collapse, then the fraction of mass that was expelled will lead to a 
one-time transient, and the proto-neutron star will accumulate the 
inner part of the envelope, perhaps transitioning to a black hole 
if it becomes sufficiently massive (in the s 1 1 _0 1 , it would form a 
black hole with a 3.85 Mq baryonic mass). But if the weak energy 
release is caused by something else, e.g. thermonuclear combustion 
of material at the surface of the degenerate (and hydrostatically- 
stable) core, the outcome is a massive eruption in a massive star. In 
this case, the luminosity eventually levels off, as the photosphere 
receding in mass through the ejecta layers eventually reaches the 
inner envelope layer that failed to eject. This second category of 
objects would be called SN impostors[f] 

Unless the event occurs nearby and allows the inspection of 
the transient site, one does not know whether a star survived or not. 
Hence, such events are ambiguous and may be interpreted either as 
a one-time non-repeatable explosion (a SN) or a potential recidivist 
like 77 Car. These are impostors in the sense that they would not 
result from the collapse of a degenerate core, but they do share that 
important property with core-collapse SNe that, in the present con- 
text, they would result from shock-heating of the envelope. Hence, 
a breakout signal would systematically herald the forthcoming ejec- 
tion, excluding its origin as a steady-state super-Eddington wind ac- 
celerated from a hydrostatic base. The longer duration of ~ 1 d and 
the softer distribution (UV-peaked) of the breakout signal should 
allow the detection of such events in large-scale sky surveys, pro- 
viding a clear discriminant between an explosive or a steady-state- 
wind origin for such mass ejections. 



3 In this case, the stellar core would eventually collapse, potentially eject- 
ing the residual envelope, which is now of much lower mass and quasi 
hydrogen-deficient. Depending on the time delay between the two events, 
this could also produce a strong interaction with the previous ejected shell, 
yielding an interacting SN. 



4.1.3 Edc-p < ^binding 

For the smallest energy deposition value in our sequence (model 
sll_0; -Bdop =5xl0 48 erg), only a tiny fraction of the envelope 
(~0.13Mq) is ejected, concomitanly with shock breakout, while 
the bulk of the envelope remains attached to the star. The en- 
ergy deposition is used to lift the envelope out of the gravita- 
tional potential well, but falling short from unbinding it. This re- 
sults in a "bloated" star whose radius has increased by about a fac- 
tor of two in 50 days after breakout, while the surface temperature 
has hardly changed. Consequently, the corresponding star, with its 
puffed-up envelope, brightens from a luminosity of ~10 5 Lq up to 
~1.37x 10 Lq, which remains well below its Eddington luminos- 
ity of 4.07 x 10 5 L Q . 

Over the two-year time-scale covered by the simulation, the 
photosphere retains a similar radius and temperature. Since the en- 
ergy left over by the shock is stored at large optical depth, and since 
it was insufficient to communicate any kinetic energy to the enve- 
lope, the transport regime is no longer of dynamic diffusion (as 
in the high energy-deposition cases) but instead of static diffusion, 
with characteristic time-scales of ~100yr for this massive star en- 
velope (Fig. |4j- Albeit weak, such an energy deposition causes a 
significant translation of the star in the HR diagram and may be 
linked, in some cases, to the photometric and spectroscopic vari- 
ability of H-rich massive stars. The central element that enables 
such low-energy perturbations to have any effect is the very low 
binding energy of massive star envelopes, in particular at the low- 
mass end. 



4.2 Effect of varying the Energy-Deposition depth 

We have argued in the introduction and in fj2]that, physically, a de- 
position of energy from deep regions in the progenitor star is more 
likely since this is where most of the energy is stored or can be 
made available. This can occur through gravitational collapse of 
the degenerate core from its original radius of a few 1000 km to 
~10km (the standard initiation of a core-collapse SN explosion), 
or from violent shell burning above the core (our preferred trigger 
for transient phenomena in massive stars). In contrast, the envelope 
layers well above the core, characterized by lower densities and 
temperatures, tend to merely react to energy changes occurring at 
depth. It is, however, interesting to investigate the dependency of 
our results on the adopted site of energy deposition. We have thus 
performed two sequences for Ed cp = 10 49 erg (model sll_l) and 
Edcp = 4 x 10 49 erg (model sll_3), with a deposition at 1.8 (as 
before), 3, 7, and 9 Mq Lagrangian mass coordinate. These cor- 
respond to locations in the helium or in the hydrogen shell. The 
deposition time in all cases is 10 s. We show the resulting light 
curves for both sequences in Fig. [10] with the corresponding key 
parameters given in Table|2] 

For the s 1 1_1 model sequence (black curve, left panel of 
Fig. ED, the energy deposition is on the order of the binding en- 
ergy of the entire envelope above the 1.8 Mq mass cut. As we 
discussed before, this led to a SN-impostor-like event, with only 
a fraction of the total envelope ejected to infinity. The envelope 
outside of a 3 Mq mass cut has, however, only a 10 47 erg binding 
energy so, as we move outward the site of energy deposition (from 
3, to 7, and 9Mq Lagrangian mass coordinate), we switch to a 
successful explosion with sizable kinetic energy and rapid light- 
curve evolution. As given in Table [2] the models with a mass cut 
at and beyond 3 Mq lead to full ejection of the layers exterior 
to the energy-deposition site, with representative ejecta expansion 
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Table 2. Same as Table [T] but now showing the results for the sequence of models with differing energy deposition depth in 
the 1 1 Mq progenitor star (see i|4.2l for discussion). 
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Figure 10. Same as Fig.[8]for the 1 1 M0inodel, but now showing the emergent bolometric-light evolution resulting for cases in which the energy is deposited 
at a Lagrangian mass of 1.8 (black), 3 (blue), 7 (green), and 9Mq (red). The left (right) panel corresponds to an energy deposition of 1 (4) X 10 49 erg. In all 
cases, we choose the time origin as the time of shock breakout. 



rates of 400-800 km s _1 and kinetic energy of 10 49 erg. Interest- 
ingly, the phase of high-brightness lasts from 90 up to 180 days 
despite the low amount of ejected material, i.e. between 1.6 and 
7.6Mq. The long duration of these events, even for low ejected 
masses, is caused by two factors. First, the expansion is relatively 
slow so that radiative diffusion is modestly facilitated by the ejecta 
density reduction with time (we are closer to a static regime of ra- 
diative diffusion). Second, the energy losses through expansion are 
minimal since we start from a very extended (i.e. loosely-bound) 
stellar envelope. Note here that in the more compact (i.e. tightly- 
bound) progenitor configurations for Type I SNe (Wolf-Rayet stars 
or white dwarfs), the cooling through expansion is tremendous and 
the events owe their large brightness at a few weeks after explosion 
entirely to the heating caused by unstable isotopes with week-long 
half-lifes. A corollary is that weak-energy ejections with little 56 Ni 
can only reach a large luminosity if they arise from a loosely-bound 
initial configuration, e.g. a big star. 

For the si 1_3 model sequence (light curve shown in the right 
panel of Fig. 1 10b. the energy deposition is greater than the binding 
energy of the envelope above the 1.8 M Q mass cut so that no matter 



what the position of the mass cut is outside of 1.8 Mq, we obtain 
full ejection of the shocked layers, with representative velocities in 
the range 500-1600 km s" 1 . 

To conclude, we find that moving the energy-deposition mass 
cut Mcut outwards for a given energy-deposition magnitude i?d ep 
produces a similar effect to increasing Sdop at a given M cut . Quan- 
titatively, this is modulated by the binding energy of the envelope 
exterior to M cu t. We observe in particular that with such RSG pro- 
genitors it is possible to obtain light curves with a very extended 
plateau even for very small ejecta mass. Note that this remains 
an experiment though, since we believe the energy deposition in 
a massive star must occur in the vicinity of the core, hence interior 
to the Helium shell of pre-SN massive stars. 

4.3 Effect of varying the Energy-Deposition duration 

A fundamental element of our study is that energy deposition must 
occur on a short time-scale to ensure the formation of a shock, 
rather than on a long time-scale which would allow radiative dif- 
fusion and convection to carry outward the extra energy. In this 
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Figure 11. Same as Fig. [8] but now showing the bolometric light evolution 
resulting from models in which the energy is deposited at a fixed rate over 
10 s (black), 1 hi' (blue), 1 d (light green), 1 week (green), and 1 month (red). 
The site of energy deposition is at Lagrangian-mass coordinate 3 Mq and 
the energy deposition is 4x 10 49 erg in all cases. The time origin is the time 
of shock breakout. 

section, we present results with V1D for a model with Edcp = 
4 x 10 49 erg and M cut = 3 Mq, but an energy-deposition time of 
10 seconds, one hour, one day, one week, and one month. Since 
the resulting ejecta and bolometric properties are similar to those 
of the model si l_3_w3 (see TableO, we do not tabulate the results 
for this sequence of models. However, we show the synthetic light 
curves m Fig.d] In all cases, a phase of shock breakout takes place, 
followed by a bright plateau that persists over a similar duration. 
No qualitative difference is visible, and quantitative differences are 
small. 

Hence, shock formation persists even for very long deposition 
times, despite the large range covered. In other words, a deposi- 
tion of energy that takes place over a time-scale much longer than 
the expected shock revival time in core-collapse SN explosion (i.e. 
one second) still leads to shock formation. We are aware that the 
neglect of feedback effects makes this exploration somewhat artifi- 
cial. In the context of a nuclear flash, expansion could, for example, 
endanger the continuation of nuclear burning. 

For energy-deposition time-scales much longer than a month, 
a diffusion wave, rather than a shock, would form, together per- 
haps with convection. Such a regime, not of interest in the present 
study, would have to be followed in 2D or 3D (see, for example, 
iMocak etal . 2008, 2009) to properly model the interplay between 
these two means of energy transport and the resulting effects on the 
envelope structure. 

4.4 Dependency on Progenitor Mass 

The choice of the 1 1 Mq model for the investigations presented so 
far was motivated by the lower density of the progenitor envelope. 
This produces larger Courant-time increments for the hydrodynam- 
ics and easier computing. However, as evidenced by the variation 
in envelope binding energy for 10-40 Mq massive stars (recall that 
these model masses correspond to those on the main sequence; the 
corresponding masses at the onset of collapse can be seen from 
Fig.|2j and by the dependencies presented in the previous sections, 
a given energy deposition will obviously produce a very different 
outcome for different progenitor masses. In this section, we thus 
explore the radiative and dynamical properties of ejecta produced 
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Figure 12. Same as Fig. [8] but now showing the bolometric light evolu- 
tion for different progenitor-mass models: 1 1 (black), 15 (blue), 20 (green), 
and 25 Mq (red; these masses correspond to main-sequence masses). The 
corresponding total progenitor masses at the time of collapse for this set is 
10.6, 12.6, 14.7, and 12.5 Mq, respectively. In all models, we deposit the 
same energy of 10 50 erg, and at the same physical location in the envelope, 
corresponding to the base of the helium shell. In the same order, this corre- 
sponds to the Lagrangian-mass coordinate of 1.8, 3.1, 5.0, and 7.2 Mq. The 
time origin is the time of shock breakout, which corresponds to a range of 
delays since energy-deposition in the four models. 

by the deposition of 10 50 erg at the base of the helium shell in 11, 
15, 20, and 25 Mq massive-star progenitors. In the same order, this 
corresponds to mass cuts at 1.8, 3.1, 5.0, and 7.2 Mq in the enve- 
lope, to envelope binding energies (we quote absolute values) of 
9.9xl0 48 erg, 4.6xl0 49 erg, 6.6xl0 49 erg, and 1.4 x 10 50 erg, and 
to total progenitor masses at the time of collapse of 10.6, 12.6, 14.7, 
and 12.5 Mq. We show the resulting light curves in Fie.[72l 

As before with the exploration of outcomes for varying 
energy-deposition magnitudes ( Ej4.lt or sites ( i]4.2l , adopting dif- 
ferent progenitor stars and depositing the energy in regions that 
are more or less bound leads to a similar modulation in light-curve 
properties. With the adopted value i5dc P =10°° erg, we obtain light 
curves reminiscent of low-luminosity Type II-P SNe if it is larger 
than the binding energy of the overlying envelope (11 and 15 Mq 
models), and long-lived and fainter events reminiscent of SN im- 
postors or low-energy transients otherwise (20 and 25 Mq mod- 
els). 

Note that if, using the same progenitors (11, 15, 20, and 
25 Mq main-sequence mass models), we deposit a given energy 
at a radius where the Lagrangian mass is 2Mq lower than the 
total progenitor mass (same shell mass underneath the progeni- 
tor surface), we obtain essentially the same light curves and ejecta 
properties for all models. This degeneracy results from the generic 
~10 47 erg binding energy of the corresponding (outer) H-rich shell 
in these different stars. 

To summarize, and perhaps paradoxically, the consideration 
of these additional progenitors (with different envelope struc- 
ture/binding energy) does not add any new perspective to our re- 
sults so far since the diversity of outcomes was fully revealed ear- 
lier on with the 1 1 Mq model through the variation of the values 
of E dep and M cut . 

Our sample of progenitor models does not include more mas- 
sive stars. Such models were not added to retain a homogeneous 
sample, focus on conceptual aspects (i.e. shock heating of a rep- 
resentative gravitationally-bound stellar envelope) and leave aside 
the specificities of the numerous different categories of stars un- 
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dergoing eruptions. One such category is that of blue-supergiant 
stars in the phase of hydrogen/helium core/shell burning. These 
objects are found at moderate masses, covered by our sample, 
but also extend to much larger masses. For example, the stars as- 
sociated with events as diverse as SN1987A and the erupting ?y 
Car are BSGs, with associated masses in the range 15-20 M0 and 
~1OOM0. However, compared to our sample of RSG stars, these 
objects are all more tightly bound, as reflected by their smaller sur- 
face radii of at most a hundred rather than a thousand times that 
of the sun. Hence, the various regimes highlighted above would 
apply with all energies scaled upwards. The duration of the tran- 
sient would depend on the amount of material ejected, reflecting 
the dependencies discussed in i]4.2l Because of their larger bind- 
ing energy, very low-energy transient resulting from the mecha- 
nism presented here would be less likely to occur in such objects. 
For the most massive objects, the proximity to the Eddington limit 
will considerably reduce the binding energy of the envelope and 
will make it more prone to eruptions if perturbed by a deep-seated 
energy source. To provide more quantitative results and compari- 
son with observations, we will present in a forthcoming study the 
gas and radiation properties we predict for very massive post main- 
sequence stars subject to shock heating. 



5 SAMPLE OF NON-LTE SYNTHETIC SPECTRA AT 5 
DAYS AFTER SHOCK BREAKOUT 

In this section, we provide synthetic spectra that illustrate the 
spectral-energy distribution of the various models discussed above. 
We choose the representative time of five days after shock breakout 
since it is a typical discovery time for SNe, and perhaps as well 
for optical transients etc. This also offers a unique reference time 
for comparison in all models. We focus on the models presented in 
^4. 1 1 in which only the energy deposition magnitude is varied, us- 
ing the 1 1 Mq model and adopting a deposition site at the 1.8 Mg 
Lagrangian-mass coordinate. 

The computation w e perform follows t he app ro ach described 
in various contexts in iDessart & Hillierl J2006h ; iDessart et al.l 
d2008l 12009b . Using the results from our V1D computations, we 
extract the gas properties at a total optical depth of a few tens and 
at a time of about five days after shock breakout. We then use 
these characteristics of radius, density, temperature, and velocity 
(see Fig.[9j( to setup the initial conditions for the non-LTE radiative- 
transfer calculation. At such early times, the photosphere resides in 
the outer hydrogen-rich part of the progenitor envelo pe, character- 
ized b y a typical RSG surface composition (see, e.g.. IDessart et al.l 
l2008h . We adopt a density exponent of ten and a homologous ex- 
pansion. Recall that the line and continuum formation regions are 
rather confined and hence such characteristics, in practice, really 
matter locally, i.e. the density/velocity di stributions may vary dif- 
ferently in a global sense (for details, see IDessart et alJl2009D . We 
show in Fig. [T3] the synthetic spectra (i.e. the emergent flux in the 
observer's frame), assuming a distance of lOMpc and neglecting 
any reddening. 

At that reference time, and from the high to the low energy- 
deposition case, we obtain a trend towards smaller ejecta velocity 
(narrower line profiles), smaller temperature (redder peak flux dis- 
tribution), combined with a smaller radius (smaller flux). At higher 
energy deposition, the synthetic spectra are reminiscent of Type 
II-P SNe (from standard to low-luminosity ones), while for lower 
energy deposition, the spectra are qualitatively similar to those of 
cool mass-losing stars like the present day 77 Car or SN impostors, 



e.g. 1961 fcwicld | l964l) . 1997b s JVan Dvk et alj|2000l) . 2000ch 
dWagner et al.||2004[) . or >002kg jMaund et alj|2006h . We thus re- 
cover spectroscopically the various regimes identified in bolomet- 
ric light curves and ejecta properties, primarily through modula- 
tions in luminosity and line width. An interesting aspect of this set 
of simulations is that for H-rich massive stars (and as long as the 
temperatures are not too low to foster the formation of molecules 
and dust), we obtain a very uniform set of spectral features (with 
the dominance of hydrogen and metals, whose abundance is that 
of the environment in which the star was born), merely varying 
in widths, with slight changes in the flux distribution mostly asso- 
ciated with variations in ionization. For example, the LBV rj Car 
has the same surface composition as the pr ogenitor of SN 1999em 
dHillier et alj200ll ; |Pessart & Hillierl2006l) ; it also has similar line- 
profile morphol ogy and flux distrib ution as that of the interact- 
ing SN 1994W l lDessart et alj|2009h . This spectroscopic degener- 
acy/uniformity adds to the difficulty of inferring the dynamical ori- 
gin of the event. 



6 DISCUSSION 

In this section, we discuss the ramifications of our results for explo- 
sions and eruptions in massive stars. Let us stress again that we do 
not make the claim that all massive-star mass ejections, or all tran- 
sients in general, stem from a strong, sudden, and deeply-rooted 
energy release in the progenitor- star envelope. But here, given this 
working hypothesis, we explore its potential for explaining the 
properties and the diversity of observed transient phenomena as- 
sociated with massive stars. 



6.1 Binding energy considerations 

The binding energy is the most fundamental quantity controlling 
the outcome of a given energy deposition. Although this conclusion 
seems trivial, its implications are far-reaching. 

Since the binding energy of massive stars increases with pro- 
genitor mass, SN explosions arising from more massive progenitors 
will, on average, be more energetic. Indeed, the energy deposition 
at the base of the envelope to be ejected has to be a few units of the 
binding energy, one unit to unbind the envelope, and the remaining 
units to give it its kinetic energy (fast expansion) and internal en- 
ergy (bright display). In the sample of massive stars studied here, 
the envelope binding energy (right panel of Fig. [3} sets a higher en- 
ergy threshold for increasing progenitor main-sequence mass. Con- 
sequently, stars originally more massive than ~15Mq require an 
energy deposition of at least 10 51 erg to unbind the envelope while 
100 times less is sufficient in a 10 Mq progenitor. 

A corollary is that quite generally, low-energy explosions in 
gravitationally-bound objects are more likely to be associated with 
loosley-bound objects. This appears counter intuitive since one may 
argue that the ejecta kinetic energy at infinity should be anything 
from zero (Edc P — ^binding) up to arbitrary large values (set by 
i?dcp — ^binding)- This is true in theory, but in practice, the range 
of energies for events classified as a SN or as a low-energy tran- 
sient differs, and vary greatly with the progenitor binding energy. 
At present, we have no reason to believe that i5do P and -EWding are 
exactly tuned in any way. This holds very obviously in the context 
of the core-collapse SN explosion mechanism where the conditions 
for a successful explosion are set by what happens in the vicinity 
of the core and thus entirely disconnected from the global thresh- 
old set by the envelope binding energy. So, we can assume that 
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Figure 13. Non-LTE synthetic spectra for the models presented in i]4. II based on the ejecta properties computed with V1D and extracted at a representative 
time of 5 days after shock breakout (we scale the synthetic fluxes adopting a distance to the object of 10 Mpc). Notice how, despite the 2-order-of-magnitude 
range in energy deposition in this model sequence, the resulting spectral-energy distributions are qualitatively similar, differing essentially in the absolute flux 
level (the larger the energy deposition, the more luminous the event) and in line-profile width (the larger the energy deposition, the larger the ejecta kinetic 
energy). This similarity is also seen in observations of hydrogen-rich transients and results from the comparable envelope composition of the progenitor star 
and the comparable ionization state at the ejecta at the photosphere. 



Edep is uniformly distributed (the number of events having an en- 
ergy between Edop and Edcp + AE scales linearly with the energy 
bin AE but is independent of £dep)> with values above and below 
the envelope-binding energy. Now, mass ejection will only occur 
if i?dep > ^binding- Consider an object with a binding energy of 
10 51 erg, and focus on the uniform distribution of events with an en- 
ergy deposition between 10 51 erg and 2x 10 51 erg (we ignore those 
events that lead to no eruption and limit the maximum energy to 
that for a typical SN ejecta, e.g. for SN1987A). We can say that all 
ejecta with a kinetic energy in the range 10 50 erg up to 10 51 erg rep- 
resent powerful explosions and will tend to be viewed/classified as 
a SN: This represents 90% of the events in our distribution. Ejecta 
with an (asymptotic) kinetic energy less than 10 50 erg will be pro- 
duced with a 10% probability, those with less than 10 49 erg with a 
1% probability, and those with less than 10 48 erg with a 0. 1% prob- 
ability. This does not mean that a low-energy transient cannot take 
place in a highly bound object, but it strongly suggests that unless 
there is a mechanism that biases the energy deposition to match the 
binding energy (all SN explosions strongly support the contrary!), 
low-energy ejecta are statistically more likely to occur in loosely- 
bound objects. Another way to consider this is that to produce a 
low-energy transient in a highly-bound object, the excess energy to 
deposit is so small compared to the binding energy of the envelope 
that the required tuning to make this low-energy ejection can hardly 



be met. Hence, low-energy transients associated with stars should 
stem from ejection of loosely-bound regions, for example the en- 
velope of low-mass massive stars or, in another context, low-mass 
shells at the surface of a white dwarf. Obviously, this is relevant 
only if any such energy deposition can take place but observations 
of transients demonstrate that it can. How it does it is a question 
left for future work. 

As more and more energy is needed to eject the envelopes of 
more massive progenitors, the extreme case of tightly bound Wolf- 
Rayet stars can only be associated with very energetic explosions. 
The fact that (nearly) all hypernovae are associated with Type Ib/c 
SNe is in this respect quite compelling. It is not so much that they 
somehow can achieve higher energy deposition; it is that they must 
(the fact that they do achieve a higher energy-deposition to explode 
is an interesting problem that concerns the explosion mechanism, 
neutrino physics, magneto-rotational effects etc). In other words, 
of all possible scenarios, only those that lead to a very large en- 
ergy deposition (well in excess of the binding energy) can yield a 
successful, and visible, explosion. This may act as a natural tuning 
ingredient to favor fas t-rotating cores in tig htly-bound progenitors 
like Wolf-Rayet stars dBurrows et alj|2007al) . The fraction of Wolf- 
Rayet stars that fail to explode following core collapse is unknown 
at present. 

This binding-energy argument may also be a very natural ex- 
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planation for the lack of detection of Type II SN explosi ons for 
objec ts with main sequence mass greater than ~20M Q dSmarttl 
l2009h . The standard neutrino mechanism may, after all, fail to de- 
liver an energy that is comparable to the binding energy of the more 
massive progenitor stars. And indeed, this may not be so surprising 
given the ten order of magnitude density contrast in the region 2- 
5 M between the 1 1 and the 25 M Q models of WHW02 (Fig.©. 
This would suggest that a lot of core-collapse events associated 
with the higher mass massive stars do not yield explosions, but 
instead f orm a black hole with no transie nt electromagne t ic dis- 
play (see lSmarttir2009h . The proposition of lKochanek et al.l d2008h 
to monitor a large number of red supergiants for failed explosion 
would nourish this option. 

It is the low binding energy (hence extended envelope struc- 
ture) of massive stars that causes the bright displays of their ejected 
envelopes, even for low energy deposition, and even in the absence 
of any heating from unstable isotopes. In general, the change from a 
modest-size progenitor star to an extended ejecta causes a dramatic 
cooling, primarily due to expansion. And indeed, in Type Ia/b/c 
SNe, this cooling is dramatic and the large SN luminosity results 
only because of the presence of unstable isotopes, whose half-life 
of days to weeks makes the SN bright on a comparable time-scale. 
In contrast, only modest radial expansion (typically by 2 orders of 
magnitude compared to 5-7 in Type I SNe) occurs after the ex- 
plosion of the already extended envelopes of H-rich massive stars. 
In such objects, large luminosities of 10 —10 L© can be achieved 
without the contribution from unstable isotopes, and even for mod- 
est ejecta kinetic energies. 

Hence, in general, low-energy (stellar) explosions/eruptions 
should occur from regions that are weakly gravitationally-bound. 
Given the very low envelope binding energy of low-mass mas- 
sive stars and hence the very low energy threshold for pro- 
ducing explosion/eruption, we argue that these should be prime 
candidates for transient events in the Universe, perhaps even 
in connection to interacting SNe. The stellar initial-mass func- 
tion is also biased in favor of such low-mass massive stars, 
which are therefore not rare (about one for every 20 solar-type 
stars) . This argumentation likely applies to SN 20 08S I Pri eto etakl 
120081: ISmith et all 12009 j; IBotticella etafl l2009h . M85-OT2006- 
1 dKulkarni et alj|2007l: IPastorello et al.ll2007l) . or the enigmatic 
SN2008ha dFolev et al.l2009l : IValenti et alj2009h . As argued above, 
the association of a low-energy low-luminosity ejecta with a Wolf- 
Rayet progenitor, believed at present to always be a highly-bound 
object, seems unlikely (although not excluded in theory) in the 
context where the energy deposition occurs deep within the pro- 
genitor star. In that context, it would require the presence of a 
loosely-bound shell at the surface of the progenitor star and a pretty 
exquisite tuning to allow the energy deposition occurring at depth 
to be no more than a few 1.001 times the envelope binding energy. 
Some of these postulates may not apply (e.g. could the energy de- 
position occur just below the surface of the WR star?), but it is un- 
clear at present which is faulty. Specific radiation-hydrodynamical 
simulations based on suitable progenitors are needed to address this 
problem. 



6.2 Explosion/eruption associated with transient phenomena 

A growing sample of transient objects possess intermediate lu- 
minosities between SNe and novae. Their origin is debated. Ow- 
ing to their low energy, the community has ar gued for either 
a weak-energy core-collaps e SN explosion (e.g. IBotticella et al.l 
120091 : IPastorello et alj[200a) . or for a super-Eddington wind in a 



SN impostor (e.g. lMaund et al.ll2006r . [Smifh et alj|2009bl) . The am- 
biguity stems from the common energy-ground in which events as- 
sociated with these two scenarios fall, the <10 50 erg kinetic en- 
ergy of the homunculus associated with 77 Car rivaling that of low- 
luminosity Type II-P SN ejecta. Reversing the argument, how many 
of these low-luminosity Type II-P SNe are in fact impostors, boast- 
ing ejecta kinetic energy that are well below that of LBVs? The 
"SN" status in this energy regime seems to be often based on faith 
rather than on unambiguous observational or theoretical evidence. 

In the shock-heated solutions we presented in EI4.1I the ejecta 
are systematically turned into a radiation-dominated plasma. At 
shock breakout, they are essentially unbound (their total energy 
Sgrav + Skin + £int is positive) and possess a large amount of 
internal energy. In those models with low-expansion velocity and 
super-Eddington luminosities of 10 6 — 10 7 Lq, the photosphere has 
properties that are comparable to those of erupting LBVs, yet their 
super-Eddington luminosity plays no role in driving the material, 
which was already unbound after shock passage. From an inspec- 
tion of the long-term light curve, it seems difficult to distinguish 
between the two since the time-scale over which the radiative dis- 
play evolves is very long and related to radiative diffusion: Apart 
from the breakout signature, it bears no imprint of the explosive 
origin, i.e. the sudden and deeply-rooted release of energy. 

The recent observation of fast-moving material ahead of 77 
Car's homunculus points very strongly to ward a shock -heated en- 
velope as the origin of the giant eruption dSmithll2008l) . The invo- 
cation of a super-Eddington wind in the low-mass progenitor star 
associated with SN2008S may be supporting the same hypothesis 
since low-mass massive stars are well below their Eddington limit, 
while a shock could easily take them beyond that limit. An inter- 
esting and generic component of all the shock-heating scenarios 
that we explored in this paper is the associated shock-breakout sig- 
nal, whose detection could therefore disentangle situations associ- 
ated with a quasi-steady wind from a hydrostatically-bound super- 
Eddington atmosphere on the one hand from those associated with 
a slowly-moving shock-heated ejecta on the other hand. 



6.3 Variability in Massive stars 

While the binding energy of massive star envelopes increases with 
progenitor mass, the mean binding energy is low for all objects that 
have retained a hydrogen envelope. In particular, and in the models 
of WHW02, that outer H-rich shell has a very uniform and very low 
binding energy of 10 47 erg. Owing to mass loss, the mass of that 
shell decreases with time for higher-mass massive stars, but it al- 
ways remains loosely bound. In our simulations, we found that even 
for a very modest energy deposition, well below the binding-energy 
value, important structural changes of the envelope occurred. 

We suspect that hydrodyna mic-fluid instab i lities taking 



place deep in the stellar interior i Bazan & Arnettl 19941 
lAsida & Arnettl l200d: lArnett et al.l [20051: iMeakin & Arnettl 



1998 



2006) 



may provide the energy seed for at least a fraction of the variabil- 
ity observed in massive stars. This has not been fully realized so 
far in part because studies of stellar variability tend to isolate the 
surface layers and ignore the interior, while studies of stellar interi- 
ors reduce all the complicated physics occurring near the surface to 
imposed and simplified condit ions at the grid outer boundary (see, 
however, ICantiello et al ,l2009h . We expect that, owing to their huge 
surface radii and therefore low binding energy, supergiant and hy- 
pergiant stars should be very sensitive to such perturbations from 
the stellar interior, and indeed obs ervations of such st ars reveal a 
rich and violent mass-loss history dSmith et al.|[2009bh . Similarly, 
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objects that through the course of their quasi-steady evolution hap- 
pen to lie close to the Eddington limit and/or reach critical rota- 
tion are naturally more exposed to potential time-dependent deeply- 
rooted energy leaks. 

This does not exclude the role of other mechanisms for mass 
ejections, for example, pulsations, rotation, radiation-driven mass 
loss etc., although it is at present not clear how these processes 
can produce the extreme properties required in the context of high- 
luminosity interacting SNe or giant eruptions of LB Vs. 



7 CONCLUSIONS 

In this paper, we have presented one-dimensional one-group (gray) 
two-temperature radiation-hydrodynamics simulations of pre-SN 
massive-star envelopes subject to a sudden release of energy above 
their degenerate core. Although at the high energy end, the likely 
phenomenon is core collapse, we more generally have in mind the 
thermonuclear incineration of intermediate-mass elements present 
in shells above the core. The motivation for this study is: 1) The 
existence of interacting SNe (which must eject a shell at least 
once before core-collapse, but more likely eject multiple shells, by 
some means not elucidated so far); 2) the identification of fast- 
moving material exterior to 77 Car's homunculus (which cannot 
stem from radiation-driving in a wind; |Smifhil2008l) ; 3) the broad 
range of energies inferred for Type II-P SNe, overlapping at the 
low-energy end with high-energy transients and massive-star erup- 
tions. The bulk of our results stem from work on the 1 1 Mq model 
of WHW02, although tests with higher-mass progenitors yield the 
same qualitative conclusions and regimes, merely shifted quanti- 
tatively. This work is an exploration of the outcome of a strong, 
sudden, and deeply-rooted energy deposition taking place at the 
base of a massive-star envelope. We are not proposing this scenario 
holds for all massive-star mass ejections, but we investigate what 
results when this circumstance obtains. There is no doubt it does, 
but how frequently and how robustly is left for future study. 

Although this result is not new, we find that the fundamental 
quantity controling the outcome of a strong, sudden, and deeply- 
rooted energy deposition is the binding energy of the stellar en- 
velope, which increases considerably with progenitor mass in the 
range 10-40 Mq. What is new, however, is our study of the long- 
term evolution of the gas and radiative properties that result from 
configurations where the energy deposited is greater, on the or- 
der of, or less than the envelope binding energy. We identify three 
regimes, with a continuous progression from 1) SN explosions at 
high energy (iSdcp > ^binding), with a complete envelope ejec- 
tion associated with a 100-day long high-plateau luminosity; 2) 
SN impostors at the intermediate energy range (iJdcp ~ ^binding), 
with a partial envelope ejection, and a more modest but longer-lived 
plateau luminosity; and 3) bloated/variable stars at the low-energy 

end (iJdcp < -^binding), w ' tn 'i tl ' e or no mass ejection but with a 
residual puffed-up envelope of weakly-modified properties. What 
conditions the results is not the magnitude of the energy deposition 
itself but how it compares with the binding energy of the envelope 
exterior to the site of deposition. Hence, to achieve the same re- 
sult requires more energy in a more massive progenitor star. These 
properties are summarized in Fig. 1141 

In all simulations presented, a shock forms systematically with 
a strength that depends on the energy deposited. It crosses the enve- 
lope in ~ 1 to ~50 days, hence communicating its energy to the en- 
tire progenitor envelope quasi-instantaneously, i.e. as opposed to, 
e.g., a diffusion time-scale for energy transport of 10 4 years or more 
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Figure 14. Correlation between various energies (normalized to the corre- 
sponding adopted energy deposition) and the energy deposition (normal- 
ized to the corresponding envelope binding energy). In black, we show the 
normalized asymptotic energy of the ejecta/outnow, and in blue (red) the 
normalized envelope/ejecta internal (kinetic) energy at the time of shock 
breakout. Only the models presented in Table 1 are shown here. The black 
curve illustrates the three regimes we have presented here, with SN explo- 
sions, SN impostors, and variable/perturbed stars as the energy deposition 
varies from being much larger, comparable, and smaller than the envelope 
binding energy. Notice how the internal energy always dominates over the 
kinetic energy at the time of shock breakout in the simulations performed 
here. 



at the corresponding depth. This shock eventually emerges at the 
progenitor surface with a breakout signal that varies from a dura- 
tion of an hour up to a few days (modulated here by the shock speed 
rather than by the atmospheric-scale height), with a flux peaking 
at longer wavelengths for weaker shock strengths. This breakout 
signal is the (and may be the only) unambiguous evidence that the 
subsequent "ejection" was triggered by shock-heating, and thus has 
an explosive origin. 

At shock breakout, the luminosity reaches a peak, then fades, 
before stabilizing/rising again forming an extended plateau. This 
plateau phase corresponds to the recombination epoch of the 
ejected mass, the internal energy decreasing primarily through ex- 
pansion and little through radiation. It is the large stored internal 
energy (in the form of radiation) that keeps the ejecta optically 
thick and luminous (radioactive decay is neglected in our work). 
We find a continuum of light curves, faster-expanding ejecta be- 
ing more luminous both at breakout (stronger shock) and during 
the plateau phase (Fig.ll5t. The models presented in ij4. ll corrob- 
orate the correlation between plat eau luminosity and mid -plateau 
photosp heric velocity identified bvlHamuv & Pinto] |2002h . and re- 
fined bv lNugent et al.U2006l) ; |Poznanski et al.l ( l2009h . We also find 
that the plateau duration is anti-correlated with energy deposition 
(Fig,|16t, At larger energy, faster expansion leads to faster cooling 
and recombination so that the ejecta photosphere recedes faster in 
mass/radius after reaching its peak earlier. For small energy vari- 
ations in this regime, interplay between kinetic and internal en- 
ergy (which are comparable at breakout) yield a plateau duration 
that is ~100d, which is on the order of Type II-P plateau lengths. 
For lower energy deposition, we switch slowly from a regime 
of dynamic diffusion to that of quasi-static diffusion. The more 
slowly-expanding ejecta, characterized by a slowly-decreasing op- 
tical depth with time, gives the bolometric luminosity a modest 
peak value and a slow evolution, with plateau durations of up to 
1-2 years. The plateau phase is thus more extended and fainter for 
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Figure 15. Correlation between the peak luminosity during the plateau 
phase (black dots) and the mass-weighted average ejecta velocity. We also 
show the correlation for the peak luminosity at shock breakout (red dots; 
scaled down by a factor of 1000 for convenience). Note that we include 
models from Tables 1 and 2. We overplot the line L oc i> 1,6 , which fol- 
lows closely the distribution of points for I/p ea k Plateau versus (v) . Our 
radiation-hydrodynami cs simulations support the correlation identified by 
lHamuv & Pintol 1 20021) and su bsequently improved upon by iNugent et al J 
j200d) : IPoznanski et alj <2009l) . Our slope is in close agreement with that 
proposed in this last reference. Impressively, the relation holds over the en- 
tire domain explored. Note that there is no consideration of radioactive de- 
cay from unstable isotopes or departures from spherical symmetry, and only 
data points associated with the 1 1 Mq -progenitor star are used. Relaxing 
these choices would likely introduce some scatter. 
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quantity equal to the asymptotic ejecta kinetic energy divided by the peak- 
plateau luminosity brightness. 



lower energy deposition, echoing the light-curve trend going from 
SNe to SN impostors (Fig.Q}. Note that, in our simulations, these 
time-scales are always 6-7 orders of magnitude larger than the time- 
scale for the energy deposition, which was chosen to be ten seconds 
in most cases. In other words, the time-scale over which the light 
curve evolves (basically that of radiative diffusion in an expanding 
medium) has no connection to the time-scale over which the energy 
was deposited in the first place. 

From our exploration and with our set of models, we find that 
explosions of varying strength can yield the broadest range of out- 



comes in low-mass massive stars because they are characterized 
by a very low envelope binding energy (Fig. [3j. We indeed obtain 
light curves evolving from week to year time-scales (Fig. [8]l and 
ejecta expansion rates ranging from a few tens to a few thousand 
km s -1 (Fig- ID- An explosion/eruption producing a transient re- 
quires merely 10 49 erg in such objects, a value that is so low that 
gravitational collapse of the stellar core may not be required. And 
indeed, in this mass range, stellar-evolutionary calculations reveal 
the existence of nuclear fl ashes in the last nuclear-burning stages 
dWeaver & Woosley|[l979l) , which could represent such an energy 
source. We therefore propose that low-mass massive stars are prime 
candidates for transient phenomena in the Univers e, as well as 
prime canditates for interacting SNe such as 1994W l lDessart et al.1 
l2009h . 

In our simulations, sudden energy deposition above the core 
leads to shock-heating of the entire envelope. Whenever the energy 
deposited is greater than its binding energy, the entire envelope is 
ejected, with an asymptotic kinetic energy that is commensurate 
with the energy deposited. If a subsequent energy deposition occurs 
(e.g. a nuclear flash followed by gravitational collapse, as needed 
in a fraction at least of interacting SNe), the second ejection would 
have low mass and little or no hydrogen. Depending on the time 
between the two ejections, one could get an interacting SN for a 
short delay (i.e. a few years) or two transients at the same location 
for a long delay (the first one being dim if the explosion energy 
is small and the second one being potentially very dim due to the 
low ejected mass). These scenarios are rather speculative, but they 
are warranted since at present most, if not all, transients have an 
unknown origin and are poorly characterized. 

In our simulations, and within the context of this work, we ob- 
tain small mass ejections only when depositing the energy close to 
the progenitor surface, an eventuality that seems difficult to justify 
physically in such massive stars. Such low-mass ejections would 
seem to be better suited, for example, to the surface layers of an ex- 
tended white dwarf. Observationally, low-mass ejections are likely 
to be associated with fast transients. For transients that are both fast 
and faint, a low-mass ejection in a highly- or moderately-bound 
object seems required. At the least, our simulations for (loosely- 
bound) RSG stars perturbed by a small deeply-rooted energy re- 
lease produce large mass ejections (the overlying hydrogen enve- 
lope) and long-faint transients. 

Synthetic spectra computed for a sequence of models 
with varying energy deposition reveal a continuous evolution 
from Type II-P SN-like spectra at high energy (L ~10 s Lq), 
to low-luminosity Type II-P SN spectra at intermediate en- 
ergy (L ~10 7 Lq), to SN-impostor-like spectra at low energy 
(L ~10 6 Lq), with, in the same order, narrower line profiles and 
redder/cooler spectral-energy distributions (Fig. [13). 

The results from this work should not be compromised by the 
approximations made in our radiation-hydrodynamics simulations. 
First, with one dimensionality we prevent convective transport and 
any structure formation through, e.g., Rayleigh-Taylor instabilities. 
This may alter the properties of models characterized by low expan- 
sion speeds (longer evolution time-scale); we thus plan to study this 
eventuality with 2D and 3D simulations in the future. Second, we 
deposit the energy at a large and fixed rate, independent of any feed- 
back effects. In the case of nuclear flashes, such feedback effects 
could shut-off the burning prematurely. We are aware of this arti- 
fact and will attempt in future work to develop a more physically- 
consistent approach by investigating the conditions that may lead 
to shock formation, rather than assuming a setup that systemati- 
cally leads to it. However, provided a shock forms, we think our 
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results apply. Third, progenitor models may differ on a mass-by- 
mass comparison with other groups but the general trend of increas- 
ing binding energy with main sequence mass should hold. Fourth, 
one-group transport should be accurate enough since it has been 
shown to ca pture the essentials of su ch radiation hydrodynamics 
simulations (Utrobin & Chugaill200^) - the key physics presented 
here takes place at large optical depth, under LTE conditions. 

Our finding that very modest energy perturbations can dra- 
matically affect the structure of a massive star motivates detailed 
multi-dimensional hydrodynamical investigations of massive-star 



collapse (see, e.g. 


. Bazan & Arnett) 1994 


19981: lAsida & Arnettl 


l200d lArnett et al 


120051: Meakin & Arnett 


2006[). As shown in 



these hydrodynamical simulations, and more generally, we sur- 
mise that the quasi-steady state approach of standard stellar- 
evolutionary codes (which keep an eye on the longer-term evolu- 
tion of stars) may be missing important ingredients for our under- 
standing of massive-star evolution. This has relevance for under- 
standing mass loss and in particular massive-star eruptions, stellar 
variability/stability, and interacting SNe. Such pre-SN mass ejec- 
tions would also modify, and perhaps sizably, the envelope mass 
and structure, thereby affecting the conditions surrounding core 
collapse and explosion. 

Ongoing and forthcoming surveys/missions like Pan- 
STARRS, Sky Mapper, the Palomar Transient Factory, GALEX, 
or the Large Synoptic Survey Telescope will better reveal the 
diversity of explosions, eruptions, and more generally transient 
phenomena in the Universe. We surmise that surveys up to now 
have missed a large number of low-energy long-lived transients, 
such as low-luminosity Type II-P SNe (objects even less energetic 
than SN1999br) and SN impostors. It is somewhat surprising that 
we have not yet detected Type II SNe with Plateau durations well 
in excess of 100 days. Moreover, for the shock-heating solutions 
presented here, a breakout signal systematically takes place. At 
SN-like energies, the signal may be too short to be resolved 
dGezari et ai]l2008l) . but for lower-energy transients, the reduced 
shock speed and strength would lengthen the breakout duration 
up to about a day, and move the peak of the spectral-energy 
distribution from ~100A to the 1000-3000A range. Hence, the 
breakout signal should be more easily detectable in such transients, 
allowing to distinguish between an explosive event and, e.g., a 
super-Eddington wind. 
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